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We study the quantum dynamics of a two-level system driven by a pulse that starts near-resonant 
for small amplitudes, yielding nonadiabatic evolution, and induces an adiabatic evolution for larger 
amplitudes. This problem is analyzed in terms of lifting of degeneracy for rising amplitudes. It is 
solved exactly for the case of linear and exponential rising. Approximate solutions are given in the 
case of power law rising. This allows us to determine approximative formulas for the lineshape of 
resonant excitation by various forms of pulses such as truncated trig-pulses. We also analyze and 
explain the various superpositions of states that can be obtained by the Half Stark Chirped Rapid 
Adiabatic Passage (Half-SCRAP) process. 
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I. INTRODUCTION 

Coherent superposition of states is a key concept of 
comtemporary quantum physics, for example in quantum 
communication and computing through entangled states 
(see e.g. Q). It is well known that various superpositions 
of states can be created in a two-level atom (states |— ) 
and |+)) travelling through a laser (or equivalently for a 
cold atom driven by a pulsed laser) or through a cavity. If 
the population initially resides in the ground state |— ), 
a laser-induced one-photon resonant process leads to a 
superposition at the end of the interaction [in the rotating 
wave approximation (RWA)] of the form 

|^) = C os [*' ffidf I-) - ze— 'sin /"*' ^-dt |+) , 
Jti 2 J t . 2 

(1) 

where Q(t) is the Rabi frequency (assumed real positive) 
proportional to the pulse envelope and the coupling, that 
is integrated between U and t / , the initial and final times 
of interaction, and lo is the optical frequency. In a one- 
photon resonant cavity, prepared in a Fock state \n+ 1), 
n > 0, the counterpart of Eq. reads 

\ip) = cos (^/n~+l ' ^Y- d ^j \-> n + 1) 

-iebi(y^+lj f ^dt\\+,n), (2) 

where \i,k) represents the atomic bare state |z) dressed 
by k photons of the cavity. For instance, a half superpo- 
sition is achieved when the Rabi frequency area is n/2. 
From the practical point of view, this process requires the 
precise control of the effectively interacting pulse shape, 
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i.e. in the case of a travelling atom, this requires a well- 
controlled and homogenous velocity and the control of 
the characteristics of the intersection between laser and 
atomic beam. As a consequence, this process is said to be 
non-robust with respect to the pulse amplitude. Creation 
of coherent superpositions of states by adiabatic passage 
is of great interest since this process requires a sufficiently 
large pulse area, but not precisely defined: hence it is ro- 
bust with respect to fluctuations of the interaction and 
also to partial knowledge of atomic and field parameters. 

The analysis of adiabatic evolution and its leading cor- 
rections is well-understood in the case of exact (one- 
photon) resonance and in the far from resonance case. 
The behavior in theses two exterme cases is qualitatively 
different. As a consequence, in the intermediate regime 
that we will characterize as quasi- resonant, there is a 
competition of different effects that to our knowledge has 
not been studied quantitatively. The goal of this article is 
to present a detailled analysis of this intermediate regime, 
whose understanding is important for practical applica- 
tions. We obtain the results combining insights from ex- 
actly solvable models with approximations obtained by 
perturbative techniques in different regimes. We test the 
quantitative validity of our results by comparison with 
direct numerical simulations. 

We can identify the different known regimes of the dy- 
namics in a two-level system considering a detuning A. 
We denote by T the characteristic time of the rising (and 
falling) of the pulse. The resonant case is defined as 

TA < 1. (3) 

We can reintepret the formulas JT} and © and extend 
them to the case of multiphoton resonant processes with 
adiabatic pulses as follows Q, H, Q . In the case of an ex- 
act n— photon resonance, the two relevant dressed states 
\—,n) and |+,0) can be considered, before the rising of 
the pulse, as exactly degenerate with respect to the dy- 
namics. The pulse rising induces a lifting of the degen- 
eracy, which leads to a splitting of the dynamics along 
the two eigenstate branches. It is important to note that 
this splitting is instantaneous at the beginning of the 
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pulse only in the case of one-photon resonance n = 1 and 
in the two-photon case n = 2 0. Higher multiphoton 
processes involve Stark shifts which modify the splitting. 
The n = 1 case induces an equal splitting along the two 
eigenstate branches. These two branches are next fol- 
lowed adiabatically by the dynamics if the pulse envelopes 
are slow enough. One can define dynamical phases which 
are the areas between the associated dressed eigenener- 
gies. When later the pulse falls, the dynamics faces the 
symmetrically inverse problem of the creation of degen- 
eracy with a recombination (instantaneous for n = 1 and 
n = 2) leading to the interference of the two branches at 
the very end of the process. The resulting final transfer 
depends on (i) the way in which the splitting (and re- 
combination) occurs (in amplitude and phase), and (ii) 
the difference of the dynamical phases. In the case of a 
complete transfer, the process has been named general- 
ized or multiphoton 7r-pulse. This process is not robust 
because the splitting (and recombination) and the differ- 
ence of the dynamical phases depend on the parameters. 
Moreover, numerics shows that it is much more sensitive 
to the detuning from the multiphoton resonance than to 
the pulse shape 0, la- 
in the opposite regime far from the resonance defined 
by the condition 



TA> 1, 



(4) 



the dynamics is at all time adiabatic in the sense that 
it follows the single dressed eigenstate whose eigenvalue 
is continuously connected to the one associated to the 
initial bare state. The nonadiabatic corrections (expo- 
nentially small for smooth pulses), that produce losses 
into the other eigenstate have been extensively studied, 
e.g. in 0. 

If the detuning is additionally time dependent (induced 
by a direct frequency chirping or by an additional off- 
resonant pulse which Stark shifts the states) and if the 
condition Q is satisfied during the rising and the falling 
of the pulse, it has been shown 0,0,0 that the topology 
of the dressed eigenenergy surfaces, as functions of the ef- 
fective time-dependent external field parameters, allows 
to determine the various possible population transfers. 
The main ingredient is a global adiabatic passage along 
one eigenstate combined with local crossings of reso- 
nances which appear as conical intersections and that can 
be precisely determined from the eigenenergy surfaces. 
This adiabatic passage results at the end of the process 
either in a complete population transfer to the excited 
state or in a complete return to the ground state. This 
process in the case of an additional Stark laser, leading 
to a complete population transfer, has been named Stark 
Chirped Rapid Adiabatic Passage (SCRAP) [l]],[ll|. 

In this paper, we study the intermediate quasi-resonant 
regimes, defined by the condition 



TA ~ 1. 



(5) 



They lead to a lifting (and creation) of a guasi- 
degeneracy. We construct formulas characterizing the 



dynamics at asymptotic times beyond the lifting of quasi- 
degeneracy, assuming adiabatic evolution along the two 
branches. 

Resonance between two quantum states and the re- 
sulting transitions are mainly understood through the 
asymptotic limit of the Landau-Zener avoided crossing 
model [T2L IT3I llil fill : A complete transition can be 
achieved by adiabatic passage beyond the avoided cross- 
ing, along the state continuously connected to the initial 
one. We will show that in the case of pulses with ampli- 
tude growing linearly in time, the problem of the lifting 
of quasi-degeneracy can be interpreted as a half Landau- 
Zener process [l6|. 

We also solve the problem of the lifting of quasi- 
degeneracy beyond the half Landau-Zener model, for 
pulses rising as power of time and as smooth exponential 
ramps. 

In the next section, we describe the model with the dif- 
ferent couplings. In Section III, we define the adiabatic 
states in the model and the conditions for adiabatic evo- 
lution along one of the adiabatic states. We show the 
dynamics when the adiabatic conditions are not satis- 
fied at early times. Section IV and V are devoted to the 
calculation of the dynamics respectively for linear and 
exponential rising coupling. In section VI and VII, we 
analyze the dynamics with perturbation theory in the 
limits of respectively large and small detuning. On the 
basis of the results of Sections IV, VI and VII, we give 
an approximative formula in Section VIII for a power law 
rising of the coupling. In Sections IX and X, we apply the 
results to obtain the lineshape of pulsed resonant excita- 
tion and to the analysis of the Half-SCRAP process. In 
Section XI, we present some conclusions and open related 
problems. 



II. THE MODEL 

We study a two- level system (states |— ) and |+)) 
driven by a near-resonant pulsed laser whose state evo- 
lution 4>(t) is given by the Schrodinger equation 



^(t)=T o H(t)0(t) 



(r) = 



B-(r) 
B + {r) 



GC 2 , (6) 



with \B-(t)\ + \B+{t)\ = 1, the scaled time r = t/T 
and the Hamiltonian in the quasi-resonant approxima- 
tion B0 



H(r) = 



-A(r) n(r) 
fi(r) A(r) 



(7) 



where we have considered the basis {|— ) , |+)}. 

We consider the Hamiltonian J7J) with a constant de- 
tuning 



A(r) = A > 



(8) 



and the following models of coupling between the initial 
Tj and a final time Tf. 
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(i) power law rising 

n(r) = ( "° n T "' T ^ n 

^ ' \ 0, T < Ti, 

with Ti = 0, Tf > 1, and falling 

q (t) = / °o (-r) n , t < r/ 
l j 1 0, r>r / , 

with Ti < — 1, 77 = 0, for an integer n > 1, 

(ii) smooth exponential rising 

Q(t) = ft e T 

with r,; — ■> — oo, 77 > 0, and falling 

Q(r) = ft e~ T 

with ^ ^ 0, 77 — ► +oo, and 
(hi) smooth Gaussian 



with 



(9) 
(10) 

(11) 
(12) 

(13) 



with Ti — > — oo and/or Tf — > +oo. We assume f2(r) > 0. 

We moreover consider for the pulse rising the initial 
conditions at time r,: B_(ji) = 1, B+(t,) = 0. 



III. ADIABATIC AND NONADIABATIC 
EVOLUTION 



-T S(t) -ry(r) 

h( T ) t s(t) 



Ha(t) = 1 

the non-adiabatic coupling 

M{t) 0(t)A(t) - n(r)A(r) 



7(t) 



and 



Since A±(t) 



(It 



A_(t) 



A 2 (r) + ft 2 (r) 



(20) 



(21) 



= R + (t)0(t) = R f (r 



(±| ^(r)) = (*±(r) 



S_(r)- 
5+(r) 

"(22) 
r)), one can in- 



terpret |A±(t)| as the population of the adiabatic states 
l ( J > ±( r ))- I n the adiabatic limit, mathematically defined 
as To — > oo, the non-adiabatic coupling 7 can be ne- 
glected and the dynamics follows the adiabatic states. 

Since Q(r) and A(r) are assumed positive, at times To 
for which the pulse goes to zero, one has 

O(ro) -> 0, (23) 

0( T o) - > and 4>a{tq) = </>(tq), hence 

A_(r ) = S_(r ), A + (t ) = B+(t ). (24) 

B. Condition of adiabatic evolution 



Adiabatic transformation 



The adiabatic states <fr±(r) are defined as the eigen- 
states of H(r), associated to the eigenvalues A±(r): 



H(t)$±(t) =A±(t)$±(t). 



(14) 



Gathering the adiabatic states in the columns of the uni- 
tary matrix R(r) = [$_(t), $+(r)] : 



with 



giving 



R(r) 



tan26>(r) 



cos 8(t) sin 6>(t) 
— sin#(r) cos#(r) 



A(r)' 



< 0{t) < tt/2, 



(15) 



(16) 



Rt(r)T H(r)R(r) 
with 



A_ 
A+ 



S(t) = VA 2 (r) + ^(r), 
we can rewrite the Schrodinger equation as 
d 

ih—4> A (T) = y\a{t)4>a(t) 



A±(t) = ±-T S(t) 
(17) 

(18) 



(19) 



In this section, we consider the rising of the pulse. 
Since one considers in the following a constant detuning 
A(r) = Aq, the dynamics is adiabatic if 



7(r) 



l7(r)| 
2T A(r) 



i.e. 



tt(r)A 



2(A 2 + ^(r)) 3/2 



< 1. 



«T , 



(25) 



(26) 



where we have defined the nonadiabatic coefficient 7(7") 
(which appears at the first order of the time-dependent 
perturbation theory). We assume here that the time 
t = Tq is a characteristic time beyond which the evo- 
lution of the system is adiabatic, i.e. without population 
transfer between adiabatic states. Our goal is to find the 
population transfer during the characteristic time before 
adiabaticity. 

For times r ~ 1 (i.e. t ~ To), Eq. I)26|) for a power law 
coupling © reads 



2T [A2 + r22] 3/2 »nrJoAo, 



(27) 



which, using £!oAo/(Ag + £!§) ~ 1, can be roughly sim- 
plified as 



2T [Aq + ng] 1/2 > n 



(28) 
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We assume in this paper that this condition (|28() is satis- 
fied for power law couplings. For large detunings defined 
here as TjAo n, the evolution is adiabatic around 
these times r ~ 1 for any Qq (and is actually adiabatic at 
anytime if one additionally excludes f2o 3> Ao, as shown 
below). For intermediate detunings defined as ToAo ~ n 
(and also for small detunings ToAo <SC n), the dynamics 
is adiabatic for r ~ 1 only when TqQq n. 

We can calculate for a power law coupling the 
scaled time tm =tu /To at which the nonadiabatic coef- 
ficient 7(t) is maximum: 



tm 



1 



2n + 1 



l/2r. 



Ao 

r» 



l/n 



(29) 



which gives the estimates 



tm = 0, for n = 1, 
™~ (^) Vn . forn>2 



and 



7(TM = Q) = i^^k- fom = l, 
l/n 



1(tm) 



(31) 



For n = 1, the nonadiabatic coefficient 7(t) decreases 



monotonically from j(tm = 0) 



i n 



to zero, with 



2 A T za 

a width of order As.. For n > 2, it is roughly bell-shaped 
and symmetric around tm- This quantity l|31l) allows to 
characterize the global adiabaticity: if 7(tm) <C 1, the 
dynamics is adiabatic at any time. This implies that 
detunings such that Ao 3> Slo induce adiabaticity at all 
times. For Ao ~ f2o, the dynamics is also adiabatic at 
all times if the detuning is large T A ^> n. (The case 
Ao ~ f2o and ToAo ~ n is not of interest here since it 
induces a nonadiabatic dynamics for t ~ To-) 

Conversely detunings such that Ao -C ^o induce a 
nonadiabatic dynamics around times tm ■ This is this last 
non trivial case, which is of interest here (accompanied 
with the condition TqQq n to have adiabaticity beyond 
t ~ 1). This case can be described in more detail as fol- 

lows: during early scaled times of order ( ^ J , the dy- 
namics is approximately adiabatic for the rising coupling 
(this initial adiabatic regime occurs only for n > 1), it 

is followed by a nonadiabatic dynamics around times tm 

/ \ l/n 

which lasts during times also of order ( Aa J (true for 

any n > 1), and an adiabatic evolution for times beyond 
(see Fig. 

For exponential and Gaussian couplings, we will con- 
sider the non-trivial case Ao <C with T)f2o S> 1, which 
allows adiabaticity [through Eq. I|26|)] for times r ~ 1. 

We will show the universality of this regime for the 
different couplings considered above. 
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FIG. 1: Adiabatic state population history |A2(r)| 2 for the 
cases of linear rising coupling n = 1 (frame a), power law 
rising with n = 2 (frame b) and n = 4 (frame c) for different 
dimensionless detuning AoTo (indicated close to each curve) 
and a fixed coupling QqTo = 100 which allows the condition 
12811 to be verified for the different cases considered. 



C. Population history 



The time evolution of the population |A + (r)| 2 = 
|(<&+(t)| 4>(t)}\ 2 in the adiabatic state |$+(t)) is shown 
on the Fig. ^ f° r power law risings. One can see the 
regions of the adiabatic and nonadiabatic evolution ac- 
cording to the analysis given below. For n — 1 the nona- 
diabatic evolution, characterized by a change of the pop- 
ulation which moves into |$+(r)), starts already for short 
times and is followed by an adiabatic evolution, charac- 
terized by a constant population in |$ + (t)) and |$_(t)). 
For n > 1 there are time intervals of adiabacity for small 
and large times. For higher n the early time interval of 
adiabaticity is longer, since the characteristic time of ris- 
ing is larger. For larger Ao the early time interval of 
adiabaticity is also longer. 
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IV. LIFTING AND CREATION OF 
QUASI-DEGENERACY BY LINEARLY RISING 
COUPLING 

The problem of lifting of quasi-dcgcncracy can be 
solved analytically in the case of linear rising of the cou- 
pling [Eq. © with n = 1], The exact solution of the 
Schrodinger equation can be found in this case, in terms 
of the parabolic cylinder functions (see Appendix^ . We 
need to calculate the dynamics at asymptotic times be- 
yond the lifting of quasi-degeneracy when the population 
of the adiabatic states is time independent. 

The initial condition A_(rj) = B_(t,) = 1, A + {ri) = 
B+(n) — leads to the amplitudes of the adiabatic states 



A±(t) ~» ^(aTfcWje*^' 1 " (32a) 
v2 



P± e 

with the transition probabilities 
1 



Ti(x±+Vd(r)) 



(32b) 



P± = 2 \aTbe t(p \ 



- (lT\/l-e-™ 2 cos^ , (33) 
^ = ar g r(l-i T J -ar g rf--z— J + -, (34) 



v2To^o 



(35) 



a = -^Vl + e-™ 2 / 2 , b = -^Vl-e-™ 2 / 2 , (36) 
v 2 v 2 

(where T denotes the Gamma-function) and the phases 
X±=Xo + arg (a =p be*) , (37) 

(38) 



'1 uj 2 \ ujW u z . 

Yn = arg r i — 1 — In — , 

v 1 2 4 / 4 4 ' 



%(t) 



\+{T)dT 



\ J^A 2 + n 2 r 2 dr. (39) 



2 ./,, 



One can interpret ^/p± exp (±zx±) as the probability am- 
plitudes of the adiabatic states from the initial bare state 
|— ) resulting from the lifting of degeneracy and the split- 
ting of the population. This splitting is accompanied by 
phases shifts ±X±. The additional phases ±rj d {r) given 
by the time integral of the adiabatic eigenvalues are thus 
the dynamical phases of the process. 

The accuracy of the asymptotics (13211 is shown in Fig. 
121 At r = 1, we observe already a precision of many 
digits both in population and phase. 

One essential result is that the adiabatic populations 
p± and the phases \± depend only on to [Eq. I|35f) ] ■ Their 
dependence is shown on Fig. |3J The phases \+ arL d X- 
go asymptotically to — 7r/2 and respectively. One can 
remark that x- is n °t very different from zero after (and 
also during) the lifting of degeneracy for any detuning. 
This trend has been numerically checked to occur for any 
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FIG. 2: History of adiabatic state population | A+ (t) | 2 [upper 
frame; numerics: full line, and analytical formula p+ 1331 1: 
dashed line] and phases arg^4±(r) [lower frame; numerics: 
full line, and analytical formulae =f(x± + 7 ?d( r )) given by Eqs. 
I|37|l and l|39|l : dashed lines] for the linear rising coupling with 
ToAo = 5 and ToQ,q — 100 (giving lo ~ 0.35). Phases are 
plotted in the interval [— 7T, +7r[, which induces artificial jumps 
that we have connected for a clearer identification. 



The reversed problem of lifting of quasi-degeneracy, 
which we called creation of quasi-degeneracy, is induced 
by a pulse falling to zero [Eq. @, with n — 1]. It leads 
to a recombination of the two adiabatic states. This has 
been calculated in Appendix lAl 



V. LIFTING AND CREATION OF 
QUASI-DEGENERACY BY EXPONENTIALLY 
RISING COUPLING 



The problem of lifting and creation of quasi-degeneracy 
can be also solved analytically in the case of exponentially 
rising and falling (|12|l coupling (see Appendix B). 
The asymptotics of the exact solution can be expressed 
in terms of the Kummer functions. In the adiabaticity re- 
gion where the population of the eigenstates is time inde- 
pendent, with the initial condition A_(r,) = B^iri) = 1, 
A+{ri) = B_|_(r,) — at t% — ► — oo, we obtain the ampli- 
tudes of the adiabatic states 



A±(t) ~» yfp± e 
with the transition probabilities 



i(«TC(T)) 



P- 



1 

l + e- 



P+ 



l + e- 



(40) 



(41) 



CO 



Scaled time x 



FIG. 3: Asymptotic adiabatic state population \A+\ 2 l|33|l 
(upper frame) and phases \- an< i X+ J37H (lower frame) as 
functions of the dimensionless quantity u) for the linear rising 
of the coupling. 

the instantaneous dimensionless pulse half-area (which is 
in fact an instantaneous Rabi frequency half-area) 

COO = y fjl^dr' =^ e\ (42) 

the dimensionless detuning 

w = T A (43) 
and the phase [given by Eq. l|B47(l ] 

C = ar g rQ+i-) +H7ln2- yln2C(r 4 ). (44) 

(In practice, is to be taken as a finite large negative 
number.) The phase £ of the amplitudes is a common 
phase for the resulting superposition of adiabatic states. 
There is no additional relative phase shift during the lift- 
ing of degeneracy. 

It is remarkable that the transition probabilities de- 
pend only on the detuning AoTo (and not on CIq). More- 
over the preceding dynamical phase is here replaced by 
a pulse area. 

The accuracy of the asymptotics (I4TJ1) is shown in Fig. 

Ill 



VI. PERTURBATION THEORY FOR LARGE 
DETUNING 

For large detuning T A , the evolution of the system 
is adiabatic at all times (we exclude f2o ^ A ). We 



FIG. 4: History of adiabatic state population p+ = jA+(r)| 2 
[upper frame; numerics: full line, and analytical formula p+ 
l|41[l : dashed line] and phases axgA±(r) [lower frame; numer- 
ics: full line, and analytical formulae £ =F £(t) given by Eqs. 
I|44|l and Q42fl: dashed lines] for the exponential rising of the 
coupling with T A = 0.4 and T tlo = 100. 



are thus interested here in small nonadiabatic correc- 
tions from the initial condition j4_(rj) = B_(rj) = 1, 
A+(ri) = B+(ri) = 0. It is well known that their main 
contribution are given by the nonsmoothness of the cou- 
pling at the beginning, characterized here by a discon- 
tinuous n th derivative 0, |2jj. In this section we give 
estimates of corrections beyond. 

It is convenient to use the adiabatic basis, in which we 
use the standard perturbation theory: 

A± (r) = A±* (t) + eA { l ] ( T ) + s 2 A^ (r) + • • • (45) 

with e = 1/ (Tq |Aq|) and the initial condition A_ (ji) — 
1, A^in) = 0, A ( ± >0 \t) = 0. The zeroth order gives 

A m {r) = e lVd(T \ Af{r) = (46) 
with the dynamical phase 

%(t) = \ I T X+[T)dT = T [ 6{T ' )dT '- m 

The first order contributions read 

£A W (t)=0 (48a) 
eA { 1\t) = ie-^W [ dT'~f(T')e 2 " ld{T '\ (48b) 
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which give at first order 



(49a) 



A+{t) « \-e-' ir > d{ ^ f dT l ^{r')e 2ir '^ T '\ (49b) 
Jn 

This shows that the phase of A- (r) is approximately 
given by only the dynamical phase. This is consistent 
with what we obtained for n = 1, where \~ jSZt was 
not very different from zero during the lifting of de- 
generacy. The transition probability at the first order 
P+ (t) = | A + (t) 1 2 at large times t — > +oo is given asymp- 
totically by 



1 



— 

4 



(50) 



For consistency of the perturbation theory, we 
should keep only terms of order e of the inte- 
gral by partial integration, using exp {2irjd{r)) — 

£ ^l + CO/Aofdr' 



exp 



We however keep here the 



full expression ij^Ujl as is usually done |Tfi| . 

We now consider the power rising coupling [Eq. JJjJ] 
with arbitrary n. Introducing the large dimensionless 
parameter 



Ao 



1/r, 



» 1 



ToAq 

and using a; = (Qq/ Aq) 1 ^ 1 t, one obtains 



with 



Jn - 2 



r 2n 



3 M " /<T ^ 1+U '^du 



dw g n (tu)e 1 ' 



defining as usually done 

ui(x) = / yl + u 2n du 
Jo 

and 



1 nx r> 



i2n 



1 + [l( W )] 

The probability to the first order reads 
P+(oo) a |J„| 2 . 



2 1 



(51) 

(52) 

(53a) 
(53b) 

(54) 

(55) 

(56) 



We are interested in the asymptotics of J„ for a n 3> 1. 
One can remark that the lower limit of integration in 
J n is xiow = 0j unlike the more standard case where 
a^iow — > — oo |l5j . This difference leads to an additional 



nonexponential contribution of the integral J, domi- 
nant for a n 3> 1, which can be calculated by n partial 
integrations: 



J n S n 



I„ 



with 



k=l 



(57) 



(58) 



and 



\lOtn) Jo 



(59) 



since g n {oo) = 0. We calculate d*~ 1 g n (0) = 0* -1 /n(O) = 
Sk,n^n\ for fc < n, giving for the nonexponential 
contribution 



1 f i 



2 la 



(60) 



The asymptotics for a n — > oo of the exponential contri- 
bution can be estimated as 



In 



with 



(-l)^e- 



J2k=o 1 F k, for even 

J2k=0 Fk 



for odd n 

(61) 



F k = (-l) fe e - Q " b " sin [^( fe +^)] e M ™ b '* cos [^(' £ +l)] (62) 



and 



/ VI - x 2n dx = — Vtt 
Jo 4n 



^ 2n) 



r( 



3w+l ' 
2n , 



(63) 



We can remark that the prefactor 7r/3 in the exponen- 
tial contribution /„ is wrong. The comparison with the 
exact result (13311 for n = 1 gives the correct prefactor 
which should be 1. This well known "ir/3— problem" is 
due to the fact that one has considered the first-order 
perturbation theory in the adiabatic basis |l5j . 

Figure [5] shows the amplitudes for n — 2 obtained with 
numerics (full lines) and with the preceding formulae 
(dotted lines). The agreement is good from ToAo ~ 5 
for the probability and from ToAo « 10 for the phase 
&rgA + . The agreement is good for arg A_ for any ui. We 
have indeed found numerically a,rgA^. « r\d for any u 
during the lifting of degeneracy. 

For a smooth lifting of degeneracy [with a coupling 
such as exponential (|llfl or Gaussian coupling (| 1 31) ] . the 
terms of type S <|60[1 are all zero. Only an exponentially 
small term contributes to the lifting of quasi-degeneracy 
for large detuning. 
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FIG. 5: Asymptotic adiabatic state population P+ = 
[upper frame; numerics (full lines), perturbation theory for 
large detuning 150|l (dotted lines), for small detuning l|7ri[l 
(squared lines) and approximate formula 1861 1. the same with 
semi-log scale plot (middle frame), and phases argvl_ and 
argA+ [lower frame; numerics (full lines), using approxima- 
tions 1491 1 (dotted lines), 17511 (squared lines) and 1851 1 as 
functions of T0A0 for the rising of the coupling with n — 2 
and Toflo = 100 (calculated at the asymptotic time r = 3Tb). 



VII. PERTURBATION THEORY FOR SMALL 
DETUNING 

For small detuning T A -C 1 it is convenient to apply 
perturbation theory with e = TqAq in the Landau- Zener 
basis as described below. The amplitudes 4>lz(t) of the 
states in the Landau-Zener basis are obtained from the 
amplitudes of the bare states via the time-independent 
unitary transformation S: 



>LZ 



Z_(r) 
Z+{r) 



with 



-71 



1 1 

-1 1 



giving the Landau-Zener Hamiltonian 
H LZ {r) =S f H(T)S 



-fi(r) -A 
™A O(r) 



associated to the Schrodinger equation 

_ d(j)LZ , 



ih- 



dr 



-(r) - T o H iz (r)0 L z(r), 



(64) 

(65) 

(66a) 
(66b) 

(67) 



FIG. 6: Same as Fig. 5 for the rising of the coupling with 
n = 3 (calculated at t — 2.5Tb). 



and the initial conditions 

z±(n) 



1 



(68) 



since one starts with A_(r,-) = S_(r,) = 1, A + (t.i) = 
B+(ri) = [see Eq. I|22|l]. The zero order of perturbation 
theory gives 



4 0) (r) 



(69) 



where £(t) is the instantaneous dimensionless pulse half- 
area 



C(r) 



To 



Q(t')(It' 



The first order contribution reads 
leading to the first order solution 



Z ± (r)^Z^(r)+eZ^(r). 



(70) 



(71) 



(72) 



For large times r — ► +00, Hlz(t) becomes diagonal and 
4>lz{t) asymptotically coincides with </>a(t), which gives 
at first order: 



A±{r) 



1 



1 + i 



T A 



dr'e ±2 < {T '' 



(73) 
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The transition probability at first order can be thus writ- 
ten as 



P+(r) 



1 



1-T A ^ dr'sin(2C(r')) 



(74a) 



For the power law rising coupling, one thus obtains for 
large times r — > +oo 



and 



1 + ^ (TK n + iL r 



(75) 



In the case of the exponential coupling (|If \ with Tj — > 
— oo, we obtain, using the variable x = XbAoexp(r), in 
the asymptotic limit r — * oo: 



j . I / „ . / " , SIM 

P+(co) » - I-ToAo / dx 



Ha) 



(81b) 



which coincides with the exact asymptotic result (|4I|I in 
the limit of small detuning. 

For the Gaussian coupling one obtains with the 

asymptotic limit taken at r = 0: 



P+(oc)« -(l-T A K n ) 



ith 



C(t) = 



2 n + l ; 



n + 1 \ "+ 



1 

TqOo 



ds sin(s n+1 ) 



( "+2 ^ 
V2(»+l)y 



and 



I 

TqQq 



[2(n + l)r^ +1 ) r (-±i, 



dscos(s n+1 ) 



2(n+l) 



/ 7r 2n+l\ 
1,2 n+1 J 



sin 



(76) 

(77) 

(78a) 
(78b) 



[ 2 ^+ 1 )r /( " +1) r(^) 



If™. 



(79a) 
(79b) 

(79c) 



For n = 1, one has in particular i^i = 0r /Xof2o/2 = 
\J tt/2uj/ToAo and 



(80) 



with w given by Eq. (|35[1 . which is recovered from the 
exact asymptotic solution (|33|l in the limit of small de- 
tuning. 

Figure shows the transition probability and the 
phases for n — 2 obtained with the preceding formu- 
las. The agreement is good until TqAq as 5. The two 
perturbative formulas for small and large detunings al- 
most match for the probability. Only a small region of 
intermediate detuning is not covered by either approxi- 
mation. 



P+(0) = -[l-G(r O )r A ] 



with 



G(x) = 



dr sin 



erf(r)) 



(82) 



(83) 



The function G{x) is shown in Fig. [7| It oscillates less 
and less for larger x, and goes very slowly to 0. We are in- 
terested in the adiabatic limit x = ToOo 3> 1, where G(x) 
depends weakly on x. This allows us to conclude that, in 
the limit of small detunings and large To^Oj the popula- 
tion transfer is quite robust with respect to the peak am- 
plitude: lim To n >i \dP + (0)/d(T o Q, )\ -C 1. Already for 
T Q n ~ 10, we have \dP+(0)/d (T Q )| < 0.05 x T A . 
Moreover we can remark that it is a bit more robust 
with respect to the detuning than the exponentially ris- 
ing case: forTo^o ~ 10, we have \dP+(0)/d (T A )| < 0.3 
(to be compared with the exponentially rising, where we 
have the constant value |dP + (0)/<ff o A o | = tt/4 w 0.79). 



VIII. APPROXIMATE FORMULA FOR 
ARBITRARY DETUNING 

The result lj7T)|) for small detuning suggests to extend 
for arbitrary n the exact results of n — 1 [Eqns lIML'l) to 
P^)l] replacing in these formulae ui by u> n : 



= - \! — K n T A 

7T 



P A C 



TqQq 



V2 



(84a) 

/ n+2 A 
\2(n+l)) 



[2(n+l)] 



n/(n+ 



!) p / 2n+l \ ' 
1 (,2(^+1) J 

(84b) 



Asymptotic analysis allows us to extend the phases as 
follows: For large ToAo, we take into account the S n 
contribution (j60|l [which gives the additional prefactor n 
in Eq. JBHJl]; for small XbA , we have used l(7H| [which 
gives the additional sine ratio in Eq. H88|) ]. We finally 
obtain 



(85) 



10 




100 



FIG. 7: Function G(x) as denned in Eq. 1831 involved in the 
population transfer for the Gaussian coupling and small 
detunings. 



with the transition probabilities 

K 1 ^ 



Pn,± 



1-4 



n,± 



1 - e" 



COS <y9„ 



I \ „ / 1 W„ 

^„ = argr -argr - - i~p 



(86) 



5. ( 87 ) 



and the phases 



Xn,— 

Xn,+ 
Xnfi 



X„,o + arg(a„ + 6„e lv "), (88) 



X«,o + narg (a„ - b n e liPn ) , 
s U 4 J 4 V 4 



(89) 
(90) 



-=VT + e-™*/ 2 , b n = -4=\/l - e"™'/ 2 
\/2 V2 



(91) 
(92) 



It is remarkable that the amplitude depends only on ui n 
(and also on the dynamical phase f) n (T) and n). This 
quantity gives the scaling of the lifting of degeneracy. 
Figure El shows the good accuracy of formula H85JI for 
any detuning for n — 2. Only a logarithmic scale shows 
that the probablities are not as precise as perturbation 
theory for large detuning. Figure displays the result 
for n — 3. The populations p 2i ± and the phase argA 2i _ 
are quite good for any detuning. The phase argA 2 .+ is 
well reproduced for TqAq < 10. 



One can remark that, as a function of the detuning, 
for T A < 5 , the phases change less fast than the pop- 
ulations. Thus the phases are more robust than the pop- 
ulations with respect to the detuning.. 



IX. APPLICATION TO THE LINESHAPE OF 
THE RESONANT EXCITATION BY STRONG 
PULSES 

In this section we apply the preceding formulation to 
calculate the transition probability after a pulse excita- 
tion of bell-like shape and of large area 



A = 2((T f ) = T J i 0(r')dr' > 1 



(93) 



to allow one to apply the preceding analysis for which a 
time where dynamics is adiabatic must occur. We con- 
sider the secant hyperbolic pulse coupling and recover the 
well-known Rosen-Zener formula. Discontinuous deriva- 
tive endings are next considered with the examples of 
truncated trig-pulses. We study the lineshape, i.e. the 
transition probability as a function of the detuning. 
We consider the initial condition B_(rj) = 1, B + (n) = 

0. 



A. General calculation 

Between the lifting and creation of quasi-degeneracy, 
the pulse is assumed to be an arbitrary smooth function 
between Tj and Tt, sufficiently far from zero to avoid in- 
termediate creation of degeneracy. The time evolution 
operator can thus be decomposed as: 



U(r f ,n) 



Ull(Tf,Ti) Ui 2 (Tf,Ti) 
= R(Tf)U c (Tf, T 2 )U a {T 2 , n)Ui(T U r 4 )R f 



(94a) 



where [see Eqs. l|Al)l. (|A30|) ] 

£/e(n,Ti) = 

and [see Eq. 

U c (Tf,T 2 ) = 



UX 2 M,Tf) [U\\(T>,T f )] 



(95a) 



(96a) 



with t 2 = —t 2 + 2tj for the power law falling and 
t' 2 = — t 2 for a smooth (exponential or Gaussian) falling, 
are associated to the evolution operators of respectively 
lifting and creation of degeneracy in the adiabatic basis, 
and 



U a {T 2 ,Ti) 



JlVd{T2)-Vd{Tl)] 







-i[ild(T 2 )-rid(ri)] 



(97) 
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with 



with the discontinuous first derivative 



Vd(r) 



To 
2 



dr 



(98) 



is associated to the adiabatic evolution between. Since 
the adiabatic states coincide with the bare states at early 
and late times, we have R'(rj) = R(r/) = 1. The ampli- 
tudes of the bare states at the end of the pulse read thus 



B-{tj) 
B+(r f ) 



Un(Tf,n), 
- [£7i2 (rf,n 



(99a) 
(99b) 



B. Secant hyperbolic pulse 



The coupling reads in this case 
f2(r) = r2osech(r) = 



H 



cosh(-r) 



with t i — ► — oo and r/ — ► +00, whose asymptotics 
rise and fall exponentially: f2(r) 2f2o exp(=Fr) for 
r — » ±00. We thus calculate B + (tj) applying Eqs. H99b|) 
and (|94bjl . and using the asymptotic result [Eqs. (|B49|) ]. 
which requires s(ti) = To J 1 J7(r)dr = To^o exp(n) 3> 
= ToAo (with ri < 0) for the lifting of quasi- 
degeneracy (see Appendix B) and TqQq exp(— t%) 
T A (with r 2 > 0) for the creation of quasi-degeneracy. 
These two conditions are satisified for Hq Aq. 
In this limit, we obtain 



B + (00) 



S_(oo) 



sin 



T fI~Sl(T)dT/2 



cosh(7rT A /2) 
sin (nT n /2) 



(100a) 
(100b) 



cosh (ttToAo/2) ' 
e 2 ^[cos (nT n /2) 

+i sin (7rT r2 /2) tanh(7rT A /2)|100c 



with £ defined in (Eq. I44fl . which allows to recover the 
well-known Rosen- Zener formula |2ll 



iB+(oo)r 



sin 2 (T n Q ir/2) 
cosh 2 (7rT Ao/2) ~ 



Pj 



RZ, 



(101) 



which is exact for any f2 and A . Note also that the 
phase of B + (oo) is also exact, but that the phase of 
B-{oo) is only approximatively valid for small TqAq. 



On 



~d~r~ 



(103) 



such that the lifting (resp. creation) of degeneracy occurs 
during the linear rising (resp. falling) of the pulse. We 
assume ri ~ 1, T2 ~ r/ — 1 in Eq. (|94b|) . We obtain the 
amplitudes of the bare states at the end of the pulse 

— cos tp cos if + sin ip (104b) 
x (iy/l - e-™ 2 + e -™ 2/2 sin^ (104c) 
B+{t p ) = -2i y /pTp~sm( X -+x+ + Vd(T p )){10U) 



. cos tp sin if — e 



-■ku> 2 /2 



sin^cos(/3 I104e) 



with p± and x± respectively defined in Eqs. 
(|3Tjl . and 



(TqAq) 2 = T A 2 
2TqOq 2 



argT 



dr~ 



T=0 



I 1 W ' 

argT( - -z— 



J33) and 



(105) 



Ji06) 



= Vd(Tp) - -0 



1 - In — 

4 



7T 

4' 



1 uj 2 
argr ( - — 



argT [l-i — 



• (107) 



Thus the lineshape |S+(r p )| is determined only by the 
dynamical phase and the first discontinuous derivative of 
the pulse. 

Figure 8 shows an example of the lineshape, accompa- 
nied with the phases, for the trig-pulse 



f2(r) = Qosinr, < r < ir, 



(108) 



giving A = 2Xof2o- We have chosen To^o = tt/2, i.e. a 
"•7r-pulse" which induces a complete population transfer 
for Ao = in this model. One can see that even for 
this 7r-pulse for which the condition of large area A S> 1 
is valid only very roughly, the numerical and analytical 
result 1104|l are very close. One can remark that as shown 
by Eq. (|104e|) . the phase of B + {t p ) is ±7r/2 and does not 
depend on the dynamical phase. 



D. Truncated pulses of power law endings 



C. Truncated pulses of linear endings 

We assume that the pulse starts and ends with linear 
rising and falling around respectiveley Ti — and Tf — 
t p > 1 : 



fi(r) 
fi(r) 



£V, 1 > t > 0, 



(102a) 



fio (7p - r) , t p - 1< t < t p (102b) 



We assume that the pulse starts and ends with power 
law rising and falling around Ti = and Tf = t p 3> 1 : 

Q(t) ~ H t", 1 > r > 0, (109a) 
n(r) ~ O (t p - t)™ , t p - 1< t < t p (109b) 

with the discontinuous n th derivative 



1 d n n 



(110) 
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FIG. 8: Lineshape |B+(r p )| 2 [upper frame; numerics (full 
lines) and using approximate formula 1104eH l. and phases 
axgA- and arg^4+ as functions of TbAo for the trig-pulse 
iTTMt and To^o = 7r/2. 



Using the generalization of the results for n > 1 of the 
preceding Section, we obtain 



B+(r p ) 



-ZtyP+P- sin (X- +X+ +%(T p ))(Hla) 



b -( t p) = P- e 



-p+e 



-i{ 2 X++Vd{T 1 



<>fllb) 



with p±, x± and rjd respectively defined in Eqs. (|86|l . 

P? |) ajid iP |l . Thus the lineshape P+ = \B+(t p )\ 2 
is determined only by the dynamical phase and the n th 
discontinuous derivative of the pulse. 

This allows to calculate approximately the lineshape 
for a trig-pulse 



Q(t) = tt sin™ r, < r < 7T, 



(112) 



where ££| = fi T™. 

cir™ It=0 u 

Figure 9 and 10 give two examples of the lineshape, 
accompanied with the phases, for the trig-pulse 



Q(t) = tt sin 2 r, < t < 7T, 



(113) 



giving „4 = T £l ir/2, respectively in "7r-pulse" and 
pulse" conditions. The numerical and analytical results 
are quite close. 



X. APPLICATION TO HALF-SCRAP 

It has been shown recently that two delayed pulsed 
lasers in an adiabatic regime can be used to yield a coher- 
ent superposition of states. This process has been named 



FIG. 9: Same as Fig. 8 for the trig-pulse fTTH) and T fi = 2 
(■7r-pulse). 
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FIG. 10: Same as Fig. 9 for T o tt = 6 (37r-pulse) 



Half SCRAP |22j. It permits to create at the end of the 
interaction a coherent superposition of states, whose am- 
plitudes do not depend on the dynamical phases and are 
thus robust with respect to the field amplitudes. One 
uses two delayed lasers: a pump one-photon resonant 
laser and an off-resonant Stark laser which allows to dy- 
namically shift the levels. The effective Hamiltonian is 
of the form Q in the basis of the dressed energies | — ; 0) 
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(state |— ) dressed by photon) and |+;— 1) (state |+) 
dressed by -1 photon), that are degenerate for the exact 
one-photon resonance (A = 0) when £1 = 0. 

We can interpret the process using the topology of the 
eigenenergy surfaces as functions of the parameters f2 
and A, combined with a local analysis of lifting (resp. 
creation) of quasi-degeneracy near the start (resp. end) 
of the process. 

The time dependence of the effective detuning A(t) = 
Ao + S(t) is only due to Stark shifts which are induced 
by the laser pulses. The process can be described by the 
diagram of the two surfaces 



A±(fi,A) =±-T V^ 2 + A 2 (114) 

which represent the eigenenergies as functions of the in- 
stantaneous effective Rabi frequency O and detuning A 
(see Fig. 1110 . The associated eigenvectors can be written 
as 



1$-) = 



cos 9 


, !*+> = 


sin# 






— sm.9 




COS0 



with 



tan26»=— , < 6 < tt/2. 
A ~ ' 



(115) 



(116) 



The surfaces display a conical intersection for O = 0, A = 
induced by the crossing of the lines corresponding 
to states | — ; 0) and |+;— 1) for O = and various A. 
We study a quasi-resonant process starting (and ending) 
close to this conical intersection. 

One essential point is that the result of the lifting (or 
creation) of quasi-degeneracy depends on the direction 
of the lifting (or creation). We can characterize two par- 
ticular directions: In the A direction with Q = 0, we 
have = 0, which means that the lifting and creation of 
degeneracy occur trivially along a unique surface: 



;0) 
-i) 



!*+>■ 



(117a) 
(117b) 



In the f2 direction with A = 0, we have 9 = 7r/4, i.e 



;0) = _(|$ + ) + !$_)) 



+;-i 



V2 



(1^ 



(118a) 
(118b) 



which implies a lifting of degeneracy occuring along the 
lower and upper surfaces with an equal sharing (1/2 if 
we consider the probabilities). In other directions, the 
lifting of degeneracy results from a competition between 
the parameters A and Q and will give a mixing of the 
two surfaces with non equal sharing in general. 

The Q direction with an arbitrary A will result in a 
lifting of quasi-degeneracy as studied in the preceding 
sections. 




LU -1 



FIG. 11: Surfaces of eigenenergies l|114|l as functions of Toil 
and To A (considered both positive). The lines in the plane 
Q = 0, connected to the respective lower and upper surfaces, 
characterize the states respectively | — ; 0) and |+; — 1). One 
example of paths involving a lifting of quasi-degeneracy is 
drawn. 



The inverse process of creation of degeneracy also de- 
pends on the direction in which the eigenstates involved 
in the dynamics arrive into the degeneracy. In the f2 
direction with A = 0, we have 

|*_) = -L(|-;0)-|+;-l» (119a) 
1*+) = i(|-;0) + |+;-l», (119b) 

To obtain a coherent superposition of states, we have 
two possibilities: 

(i) first lifting of degeneracy in the A direction with 
= [according to Eq. Q117a)) ] giving one single dressed 
state involved in the dynamics; next adiabatic following 
on this dressed state (along the lower surface) and finally 
creation of quasi-degeneracy in the O direction [according 
to Eq. (|119a|) for the case of exact resonance A = 0] ; 

(ii) first lifting of quasi-degeneracy in the Q direction 
[according to Eq. I|118afl for the case of exact resonance 
A = 0] giving two dressed states involved in the dynam- 
ics; next independent adiabatic following on these two 
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dressed states (along both the lower and upper surfaces) 
and finally creation of degeneracy in the A direction with 
ft = [according to Eq (fTTTjl ]. 

These two cases are produced by two following differ- 
ent sequences of pulses: respectively (i) first the Stark 
pulse and next the pump pulse (referred to as Stark- 
pump sequence), (ii) first the pump pulse and next the 
Stark pulse (referred to as pump-Stark sequence). (Both 
sequences require an overlapping of the two pulses to in- 
duce adiabatic following.) 

The dynamical phases coming from these two differ- 
ent sequences are not identical. For the sequence Stark- 
pump at exact resonance A = 0, we start with the lifting 
of degeneracy \4>{ti)) = |$-), followed by an adiabatic 
passage 



\c/ ) (t))=e~ iI ^ dsX - is)/h 
which leads at the final time i/ to 

l^(t / )) = -L e - i 4 /dsA -w/ ft (|-)- 



$-) 



(120) 



-)) (121) 



(where a coherent state for the photon field has been 
considered.) For the sequence pump-Stark at exact res- 
onance A = 0, we start with the lifting of degeneracy 
\<j) (ti)) = (|$+) + |$-)) /\/2. The dynamics is next char- 
acterized by an adiabatic passage along each branch: 



i 

71 



-if l dsX + (s)/h 



-iJl ds\-(s)/h 



which leads at time tf to 



(122) 



|0 (t/)) = -ie"'^ JsA+(s)/fi 
v2 



ot | + ) +c i/ t t / f (A+M-A-M) 



Thus the two sequences (with exact one-photon reso- 
nance) lead to the same probabilities 1/2 but with differ- 
ent phases. The pump-Stark sequence (|123fl leads to a co- 
herent superposition of states with a (non-robust) phase 
difference f* ds (A+(s) — A_(s)) /H, coinciding with the 
dynamical phase difference, in addition to the optical 
phase. 

If one considers non exact one-photon resonance A ^= 
0, one obtains, by lifting of quasi-degeneracy, additional 
phases and amplitudes different from 1 / V2, that will de- 
pend on the shape of the pulses according to the analysis 
of the preceding sections. Fig. El gathers populations of 
different coherent superposition of states given by half- 
SCRAP for various pulse shapes. One can see that ro- 
bustness with respect to the detuning is better for power 
law rising than for the exponential and Gaussian rising 
(better for small n and large fio)- Robustness with re- 
spect to the amplitude f2o is better for smoother rising 
(and better for the exponential rising, which is indepen- 
dent of fioi than for the Gaussian rising, which is weakly 
dependent on Qq). 
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FIG. 12: Half-SCRAP population transfers P+ = |(+|0(t/))| 2 
for power rising of the pump laser (upper frame: n = 2, mid- 
dle frame: n = 4) with respectively Tof2o = 1000, 100, 10 from 
(upper) left to right (indicated close to each curve) [numer- 
ics: full lines, dashed lines: formulas (1861 1 . and for Gaussian 
rising with respectively Tbfio = 1000, 100, 10 and exponential 
rising from (upper) left to right (lower frame) as a function 
of T A . 



XI. CONCLUSION 

In this article, we have analyzed the dynamics associ- 
ated with lifting of quasi-degencracy with a constant de- 
tuning. This is the situation encountered for one-photon 
quasi-resonance. The dynamics becomes quite different 
for n— photon quasi-resonance in effective two-level mod- 
els, which leads to an effective time-dependent detuning. 
Indeed, in this case the effective Hamiltonian is, in the 
basis {| — ; 0) , | + ; — n)}, of the form 



H(r) 



-a£ 2 ( T )-A I3£ n {r) 
P*£ n {r) a£ 2 (r) + A 



(124) 



with £{t) the field amplitude, a real (chosen positive), 
(3 = |/3|e lv , and where only the leading second order 
(Stark shifts) have been kept in the diagonal. If n = 2, all 
the terms of the matrix involving the field amplitude have 
the same order which complicates the lifting of quasi- 
degeneracy. If A = 0, we obtain the following lifting of 
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degeneracy 0,3 

|-;0) = sin0|$ + )+cos0|$_), (125a) 
\+;-n) = e lip cos6\<f> + ) - sin0 , (125b) 

with 

1/31 

tan20=^, < 6» < tt/2. (126) 

a 

If n > 2, it is known that "one can approximately com- 
pensate the Stark shift" . This can be explained and ex- 
tended as follows on the example for the rising of the 
pulse: we separate the dynamics following the different 
orders in £ of the matrix: (i) the Stark shift at early times 
shifts the diagonal elements (i.e. the dressed states in an 
effective way) for small field amplitudes without transfer- 
ring any population, and (ii) for larger field amplitudes, 
the coupling lifts the resulting quasi-degeneracy. This 
last step can be in principle analyzed by the tools pre- 
sented in this article. Thus, by adjusting A, it is possible 
to cancel almost completely or more generally partially 
the effect of the Stark shift, leading to an arbitrary (in 



probability) superposition of states. The complete can- 
cellation is for example the key that permits to orient 
molecules by adiabatic passage (in that case the Stark 
shift is of exponential order) |25|. 

The analysis of the lifting of quasi-degeneracy is more 
complicated in the case n = 2 since the Stark shifts and 
the lifting occur simultaneously, leading to a lifting of 
quasi-degeneracy with a time-dependent detuning. This 
requires the extension of the tools presented here. 

The present analysis has allowed to recover the line- 
shape for secant hyperbolic pulses (Rosen-Zener formula) 
and to determine quite precisely the lineshape for trig- 
pulses. The lineshape for Gaussian pulses is to our knowl- 
edge an open question. 
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APPENDIX A: EXACT SOLUTION FOR A 
LINEARLY RISING COUPLING 

The linearly rising (resp. falling) coupling problem 
can be solved exaclty by expressing it as a Landau- Zencr 
problem of finite duration, with a start (resp. end) at the 
avoided crossing, the so-called half Landau-Zener prob- 
lem |16|. In terms of time evolution operators, we have 
to solve 

m-^U(r,T )=T H(r)U(T,T ), U(t ,t ) = 1 (Al) 
and in the adiabatic basis 
ih-^-U A (T, t ) = H a (t)U a (t, t ), U a (to,t ) = 1 (A2) 

OT 



with 

U a (t,to) = RHt)U(t,t )R(t q ) 
which can be written as 



Ua(t,t ) 



U\\t,to) U a Ht,t ) 
-[U\ 2 (t,t )Y [U£(t,t )Y 



(A3) 



(A4) 



[property also true for U(r, To)], since the Hamiltonian 
H(t) is of trace 0. (* stands for the complex conjugate.) 
We consider the linearly rising coupling 



n T, t > 
0, t < 



(A5) 



starting at t = Tj = 0. We will denote U a+ (t, t ) [resp. 
U A - (r, to)] the evolution operator asociated to the rising 
(resp. falling) of the coupling 



Landau-Zener basis 



The amplitudes </>d(t) of the states in the Landau- 
Zener basis are obtained from the amplitudes of the bare 
states via the time-independent unitary transformation 



>LZ 



with 



S = 



Z_(r) 
Z + (t) 



= S f 0(r), 



V2 



1 1 

-1 1 



(A6) 



(A7) 



giving the Landau-Zener Hamiltonian 
H lz (t) = S+H(r)S 

_ n r -q(t) -a q 

_ 2 [ -A fi(r) 
We introduce the dimensionless time variable 

T(t) = 



that leads to the Schrodinger equation 

9(j)LZ 



ih'-^-{T) = W L z{T)^ LZ {T), 

<I>lz(t) 



OT 
kz(T) 



Z-(T) 
Z+(T) 



with the Hamiltonian 

Hlz(T) = h 

the dimensionless coupling 



—T -w 
-u +T 



T A 



\/2Toflo 

and the initial conditions at time 

T t = T{ n ) = 0. 

2. Exact solution 



(A8a) 
(A8b) 

(A9) 

(AlOa) 
(AlOb) 

(All) 

(A12) 
(A13) 



The problem is reduced to the finite Landau-Zener 
model with a start at T — Tj = 0. Using the results 
of Ref. 0| one can write the exact solution of this prob- 
lem: 



0Lz(T) = U LZ (T,T l )0 LZ (T i ) 



with 



U LZ (T,T z )=^U(T,n)S 
and the evolution operator 



(A14) 



(A15) 



U LZ (T,T t ) : 
whose matrix elements read 



U^ z {T,T t ) Ul%(T,Tj 
-[Ull(T, Tl )Y [Ul z (T, Tl )Y 



(A16) 
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Ujf z {T,Ti) 



/2tt 



1+iw 2 / 2 I ?i 



-l+iw 2 /2 



r(i- 



1;, 



i7r/4 



(r^ e - 47r/4 ) Ac V 2 (^V2 e ^/ 4 ) 



(Al7a) 



(A17b) 



where D v (z) represents the parabolic cylinder function of order v and argument z. Since in our case Tj = and using 



D v (0) 



f(¥) : 



we obtain 



ur z (T,T it 



r (1 - fWu; 2 ) 2-w 2 / 4 

r (i - 

2^ 2 / 4 e^/ 4 r (1 - \iu 2 ) 



72 



iw 2 /2 



(rV2e~ 47r / 4 ) + Au,V2 (rV2e 43 ^ 4 
(rV2e-™/ 4 ) - A^ 2 /2 (rV2e i3 ^ 4 ) 



(A18) 

(A19a) 
(A19b) 



a. Asymptotics 



We calculate the asymptotic expansion for the evolution matrix Ulz when the evolution is adiabatic, i.e. for r > 1 
and for TqFIq 3> 1 since we will consider intermediate or small detunings. We will thus determine the asymptotic 
expansion in the limit T> 1. Following Ref. [lfij . we have the leading terms of the large-argument and large-order 
asymptotic expansion 



D 



iuj 2 /2 



D 



iu 2 /2 



(TV2e~ tw/4 
(TV2e i37r/4 



with 



cos$(T) exp 
cos$(T) exp 



-u> +ir) 
3tt 2 

— —UJ +11] 



■ sini?(T) exp 



-~U - l [7]+ - 



(A20a) 
(A20b) 



(A21) 



and 



V 



1 In 

4 2 



#(T) = --0{t), 



4= (T + Vw 2 + T 2 



+ | V / ^ 2 +T 2 . 



(A22) 



(A23) 



This asymptotic expansion (|A20|I has been checked numerically to be in fact valid when either to or T is large. This 
will allow us to use it for any detuning Ao in the adiabatic region, i.e. when T> 1. Moreover this asymptotic 
expansion is already a good approximation for T ~ 3. Using 
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and 



we obtain 



with 



r(z)r(z + l/2) = ^r(2z), zeC > (A26) 

U^zi^Ti) ~* ae i7)l(T) costf(T) + &e~"' 2(T) sintf(T) (A27a) 
Uf z {T,Ti) ■** fce^^cosi^T) - ae""' 1 ^ sintf(T) (A27b) 



TftfT) = ^(r)-^(^l-ln^+argrQ-^, (A28a) 
m(T) = UT)-^(l-^)+l+^gT(l-^iw 2 Y (A28b) 



rj d {T) = f V^ 2 +T 2 dT=^-V^ 2 + T 2 + %-\n T+Vu;2 + T2 , (A28c) 
Jo 2 2 uj 

= % (r) = i / A+(r)dr (A28d) 







= -^=^1 + 6-^/2, b= -Ly/l-e-™ 2 / 2 , 
V 2 v 2 



(A28e) 



where TjdiT) corresponds to the dynamical phase associ- 
ated to the positive instantaneous eigenvalue Vw 2 + T 2 
of the Landau-Zener Hamiltonian V\lz{T)/Tl of Eq. 
(IATT1) . 



3. Time evolution operator for the adiabatic states 

The evolution operator in the basis of the adiabatic 
states can be written as 

U a+ (t, 0) = R f (r) SU lz (T(t), 0)S^R (0) , (A29) 

with R(0) = 1. Thus, 



ae ^i(T(r)) +5e ^ 2 (T(r)) 
_ oe *Ji(T(T)) + &e ir, 2 (T(r)) 



(A30a) 
(A30b) 



APPENDIX B: EXACT SOLUTION FOR 
EXPONENTIALLY RISING COUPLING 

The case of exponentially rising and falling coupling 
can be solved analytically in terms of the Kummer func- 
tions. Here we give the asymptotics of the evolution op- 
erator in the adiabaticity region where the population of 
the eigenstates is time independent. 

We consider the rising coupling 



with the initial condition at r$ — > — oo. 



(B31) 



1. Evolution operator for the bare states 
It is convenient to introduce the new variable [23l 124 
s(r) = T f n{r')dr' = T Q e T , 

J —OO 



(B32) 



that corresponds to the partial dimensionless pulse area. 
In terms of this new variable, the Schrodinger equation 
reads 



-0(a) 1 
1 6(s 



B+(s) 



= <Kt) = 



B-(r) 
B + (r) 



where 



Q(s) 



A[r( 8 )] 
Q[r(s)] 



(B33a) 
(B33b) 

(B34) 



is called the Stueckelberg variable. In our case A [r (a)] 
A , A [r (a)] - a/To, i. e. 



9(a) 



T A 



(B35) 



The differential equation for -B-(a) for an arbitrary func- 
tion (a) reads 



-2i6(a) + e(a) 2 + l B_(s), (B36) 
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which gives for the exponential coupling 



d 2 B_ 
ds 2 



2i 



T A , (T A ) 2 



S_( S ). 



(B37) 

One starts from the time Tj going to — oo. It means that 
Si = T()D,Qe Ti goes to zero (positive). The general initial 
conditions for -B-(s) are 



— (a i ) = -\B + (n)-—B-{7 i ) 

We introduce the new function 

C-{s) =B-{s)s- im ' 2 

with 

TD = T Aq, 

whose evolution equation is 



d?CL 
ds 2 



(s)+i 



dC„ 



1 



s ds + 



The general initial conditions for C_(s) read 

C-(s i )=B_( Ti )s7™/ 2 , 
dC. 



(B38a) 
(B38b) 

(B39) 
(B40) 

(B41) 

(B42a) 
(B42b) 



Considering the particular initial condition B-(ri) = 1, 
B + (Ti) = will give the components fn(r, Tj) = -B-(t) 
and U 2 i{t, Tj) = B+(t) of the evolution operator {/ (r, Tj), 
which is enough to characterize completely the evolu- 
tion operator since we have ^22 (t, Tj) = U*i(r,Ti) and 
U\ 2 (t, n) = — C^2i( t j T i) (having the Hamiltonian H(t) of 
trace 0). 

Eq. (|B41|I is a confluent hypergeometric equation, 
which for the initial conditions C_(sj) = lro / 2 ; 
dC, -( s ') = has the solution 

ds 

C_ (s) = s ; m/2 e~ ls/2 M(izu/2, iw, is), (B43) 
where M(a,b,z) represents the Kummcr function. 

2. Asymptotics 

We are interested in the limit Ao <C ^0 with ToSlo ^> 1, 
which corresponds to the asymptotics for large s 3> w: 



M(izu/2,iw,is) (is) 



-iva/2 r (izu) 

T(im/2) 



(e ls + e -* m/2 } . 

(B44) 



Using Eq. (|A26(I we obtain 











a- 




cosh ~ 


-I) 











(B45) 



This leads to 

C7ii(T,7-i) = [U 22 (T, Ti )]* = B_(r) - (B46a) 

(B46b) 



pis/2 _i_ „-(7rro+is)/2 



v /2(r+e =¥ ^) 
l/2i (t,t0 - - [U 12 {T,n)]* = B+(t) = S+(s) (B46c) 

e -(7riu+is)/2 _ gis/2 



~~> e 



V2(r 



(B46d) 



with 



i ' ( +roln2-^ln Si) (B47) 



where we have used the definition of B + (s) : B + (s) = 
2i«l + HB_( s ) and 

as s V J 



|r(l/2 + ia)| = 



a £ 



(B48) 



cosh (ita) ' 

3. Evolution operator for the adiabatic states 

The evolution operator in the basis of the adiabatic 
states U A (T,n) = R 1 '(t)C/(t, t ; )R(t 1 ), R(r. t ) = 1, has to 
be considered for r ~ 1 and Ao <C ilo which gives = 
7r/4 for the transformation R^(t), according to Eq. I|16|l . 
This leads to 



U%_(r,Ti) = -Wtfii(T,7i) - CT 2 i (T,Ti)) (B49a) 



1 



V2 

e i(£+s/2) 



vr 

1 



(B49b) 



Uf + {T,Ti) = —(Uu(T,Ti) + U 21 (T,Ti)) (B49c) 



V2 

e -7rw/2 e i(£-s/2) 
^ +e -™ 



(B49d) 



If one starts with the initial condition A_(— 00) = 1, 
A + (— 00) = 0. the amplitudes of the adiabatic states are 
A_(r) = ^(r.rO and A + (t) = C/f^n). 



APPENDIX C: FALLING COUPLING AND 
CREATION OF DEGENERACY 

Using the time reversal symmetry, we calculate in this 
appendix, from the assumed known evolution operator 
Ua+(t, Tj) characterizing the lifting of degeneracy start- 
ing at t = Tj, the evolution operator Ua— (t/,t) charac- 
terizing the creation of degeneracy. We assume a process 
with a falling coupling starting in the adiabatic region 
and ending at r = Tf. The coupling falling as a power 
law HI Up starts at t < — 1 and ends at t = Tf — 0; the ex- 
ponentially falling coupling i|12fl starts at r < and ends 
at t = Tf — > +00. The Hamiltonian associated to the 
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pulse falling H_(r) is determined from the Hamiltonian 
associated to the pulse rising H + (r) by H_(r) = H + (— r). 
The relation between the adiabatic Hamiltonians is 



H a _(t) = [H a+ (-t)]' 



(C50) 



with H J 4 + (r) and H A —(r) the Hamiltonian in the adi- 
abatic states associated respectively to the rising and 
falling of the coupling. The Schrodinger equation reads 
in this case (for any to) 

ih^-U A ^(T,T Q ) = H A _(r)[/A-(r,r ). (C51) 
or 

After complex conjugation, using r' = — r, we obtain 

U a _(-t',t ) = [U a+ (t',t \)}* (C52) 
with Tq = —To in order to satisfy U A -{tq, tq) = 

[f/A+(-TQ,^)]* = 1. ThuS 

U A -(t ,-t') = [U a .(-t',t )] 1 = [U a+ (t',-t )}\ 

(C53) 

where t stands for the transposed. Choosing tq = Tf, we 
get for any r in the adiabatic region 

U A -(tj.t) = [UA+i-^-Tf)]* (C54a) 
U A \(-T,-T f ) -[U\\(-r,-r f )\ 



U\\(-T,-T f ) [U\\(- 



-v 



)Y 

(C54b) 



with Tf = for the coupling falling as a power law i|10|) 
and Tf — > +00 for the exponentially falling coupling (|12H . 
If the coupling falling as a power law ends at r/ 7^ : 



n( T ) = /- fi o(r-r / ) B ) r<r/ 
w \ 0, r > T f 

we obtain using the same method 



(C55) 



with 

I n=(— ) / +0 °^ 5 „He— . (D3) 

Thus standard techniques of contour integration allow to 
calculate with high accuracy Im/„ for n even and Re/„ 
for n odd as follows. We have to evaluate 



-1 



(D4) 



In the upper half complex plane, g n [ui{x)\ has n poles 
(with x extended in the complex plane) 

x (k) = e i%(i+k) ^ fc = o,l,...,n-l, (D5) 

that are associated to poles 



(fe) 



\/T + u 2n du 



b x^ 



of order one in the variable ui with 



(D6a) 
(D6b) 



b n ee jf Vl-x^dx = i-V^^|4iy. (D7) 

We evaluate the function g n (ui) around each of these 
poles as 



1 x 



(fc) 



fi ( fc ) 

U) — Ulc 



using the relations around the poles 

-2n 



1 + x* 



„(fc) 



(D8) 



(D9) 



U A -(Tf,r) = [U A+ (-T + 2Tf,Tf)] t 



(C56) and 



APPENDIX D: ASYMPTOTICS FOR LARGE 
DETUNING: COMPLEMENTS 

In this appendix, we show formula (|fc> 1 1> - starting from 
the integral JSHl 



I n ee — / 92ffn(w)e l 



(Dl) 



We first notice that d™g n (ui) are odd functions for any n. 
Thus 

Im/„ = —Iml' n for n even, (D2a) 
Re/„ = ^Rel' n for n odd, (D2b) 



3/2 



This simply leads to 



«sf >M » <-i)" +1 ^ 



(*>)• 



6 / (fc) 



n+1 



Hence, d^g^\uj) having n poles of order n+1, we obtain 
I' n =2m J2 Res [dZ9»(u)e ia ^] (D12) 

lm(w)>0 

with the residue for each pole 



(D10) 



(Dll) 
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Res [dZgnHe 1 ' 



x lim 9,™ 



(D13a) 
(D13b) 



r 



giving 



lb — ± 



This leads to 



fc=0 



(D14) 



x Efe=o 1 * Im i for even n 



2tt 
3 



2 e 



with 



ft 



for odd n 
(D15) 



-)fc e -a„fc„sin[^(fc+i)] e ia„fc„cos[-I-(fe+i)]^ (D16) 



We can remark that this result is approximative due to 
the approximation (|Dlljl . It is more accurate for larger 



a n . This allows to calculate Im/„ for n even and Re/ n 
for n odd as prescribed by i|D2|) . We have additionally 
found numerically that extending I' n as follows 



I" 



2t 
3 



K-i)- 

and setting for all n 



^ElLVFu, for even 

— a n b n 



J2k=0 ^ k 



-I" 



, for odd n 
(D17) 



(D18) 



which leads to Eq. I|61(l . allows to calculate /„ with a 
good approximation for sufficiently large values of a n . 



